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Abstract 

Vector  plane-wave  superpositions  defined  by  a given  set  of  orthonormal  scalar  functions  on 
a two-  or  three-dimensional  manifold  are  treated.  We  present  a technique  for  composing 
orthonormal  electromagnetic  beams  in  complex  media  and  free  space.  New  families  of  exact 
time-harmonic  solutions  of  the  homogeneous  Maxwell  equations  in  linear  media  (isotropic, 
chiral,  anisotropic,  and  bianisotropic)  and  free  space,  which  describe  orthonormal  electro- 
magnetic beams,  are  obtained. 


1.  Introduction 

Properties  of  electromagnetic  beams  in  composite  and  complex  media  (chiral,  anisotropic,  bian- 
isotropic) are  still  investigated  insufficiently.  In  some  practical  cases,  for  example,  in  character- 
izing complex  media  by  free  space  techniques,  the  plane-wave  approximation  of  incident  beams 
proves  to  be  inadequate.  As  in  the  case  of  conic  refraction  in  biaxial  crystals,  in  investigating 
possible  physical  phenomena  caused  by  anomalous  propagation  in  homogeneous  [1,  2,  3,  4]  and 
helicoidally  nonhomogeneous  [5,  6,  7]  bianisotropic  media,  the  plane-wave  model  of  the  incident 
beam  is  also  inadequate.  The  most  commonly  used  Hermite-Gaussian  and  Laguerre-Gaussian 
beams  are  mere  the  approximate  solutions  of  Maxwell’s  equations  in  free  space. 

In  Ref.  [8],  a novel  technique  for  composing  orthonormal  electromagnetic  beams  and  some 
other  specific  exact  solutions  of  wave  equations  in  linear  media  is  suggested.  By  applying  this 
technique,  some  new  families  of  exact  solutions  of  the  homogeneous  Maxwell  equations  for 
electromagnetic  waves  in  isotropic  media  and  free  space  are  obtained  in  Refs.  [8,  9,  10].  In  this 
paper,  we  briefly  outline  the  basic  ideas  of  this  technique  and  present  its  further  applications. 


2.  Beam  Manifold,  Beam  Base,  and  Beam  State 


Let  ( un ) be  a set  of  complex  scalar  functions  un  : Bu-tCl  on  a real  manifold  Bu,  satisfying  the 
orthogonality  conditions 

{um\un)  = [ um{b)Un(b)dB  = Smn,  (1) 

JB  u 

where  dB  is  the  infinitesimal  element  of  Bu;  is  the  complex  conjugate  function  to  um,  Smn 
is  the  Kronecker  ^-function.  Let  us  consider  a plane-wave  superposition  (termed  below  the 
“beam”  for  brevity  sake) 


Wn(x)  = 


efx-K(6) 


ix-K(6) 


un(b)u(b)W{b)dB 


un(b)u{b)W{b)dB, 


(2) 
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where  B C Bu — beam  manifold — is  a subset  of  Bu  with  nonvanishing  values  of  function  W'  = 
v{b) W(6);  x and  K are  the  four-dimensional  position  and  wave  vectors.  For  electromagnetic 
waves,  W can  be  any  of  the  following  quantities:  the  electric  (magnetic)  field  vector  E (B),  the 
four-dimensional  field  tensor  F,  or  the  six-dimensional  vector  col(E,H).  Functions  K = K (b) 
and  W = W (b)  specify  the  set  of  plane  harmonic  waves  involved  in  the  beam  (beam  base), 
whereas  a complex  scalar  function  v = u(b)  specifies  the  beam  state.  The  function  Wn(x)  is  an 
exact  solution  of  the  linear  field  equations,  provided  that  the  integral  in  Eq.  (2)  exists. 

The  fields  treated  in  this  paper  are  composed  of  plane  waves  with  factorized  vector  ampli- 
tudes W (b)v(b)un(b)dB,  where  each  of  the  factors  plays  a specific  part  in  determining  the  field 
structure.  The  normalized  vector  factor  W(6)  prescribes  the  polarization  of  the  infinitesimal 
plane  wave,  whereas  the  scalar  factors  specify  the  intensity  and  the  phase  as  follows.  The  func- 
tion v = u{b)  is  used  mainly  to  obtain  a set  of  orthonormal  beams  or  to  change  the  beam  state 
[8].  In  some  special  cases,  it  simply  reduces  to  a normalizing  constant  factor.  In  the  general 
case,  the  infinitesimal  element  dB  depends  on  b,  so  that  both  un  and  dB  act  as  weight  functions 
defining  the  ratios  in  which  plane  waves  with  different  propagation  directions  and  polarization 
states  are  mixed  in  the  beam  W„.  However,  the  distinctive  features  of  Wn(x)  as  opposed  to 
other  beams  are  imprinted  completely  by  the  complex  factor  un  = un(b).  This  factorization  is 
equally  convenient  both  to  set  a unique  wave  pattern  for  each  beam  and  to  obtain  the  whole 
family  of  beams  with  a prescribed  general  property  such  as  the  orthonormality. 


3.  Time— Harmonic  Orthonormal  Beams 


There  are  four  key  elements  defining  the  properties  of  the  presented  beams:  the  set  of  complex 
scalar  functions  un  = un{b)  on  a real  manifold  Bu,  the  beam  manifold  B,  the  beam  base  functions 
K = K(6)  and  W = W(6),  and  the  beam  state  function  v = i /(b).  By  setting  these  elements 
in  various  ways,  one  can  compose  a multitude  of  normalized  and  orthonormal  beams  with  very 
interesting  properties  [8,  9j.  In  particular,  these  elements  can  be  set  [8,  9]  in  such  a way  as  to 
satisfy  the  orthonormality  condition  smn  = (Wm|Q|Wn)  = NQ6mn,  where  Nq  is  some  given 
scalar  coefficient. 

A time-harmonic  beam  of  the  form  Wn  [Eq.  (2)]  with  a two-dimensional  beam  manifold  B, 
propagating  in  a homogeneous  linear  medium,  can  be  written  as 

Wn(r,  t)  = e~iult  [ eir  *Wun(b)v(b)W(b)dB.  (3) 

Jb 

Let  us  introduce  a scalar  product 

Smn  = (Wm\Q\Wn)=  f W^(r, t)QW„(r,  t)dao,  (4) 

J<T0 

where  <jq  is  the  plane  with  unit  normal  q,  passing  through  the  point  r = 0,  Q is  some  Hermitian 
operator,  and  W}n(r,  t)  is  the  Hermitian  conjugate  of  Wm(r,  t).  To  compose  orthonormal  beams, 
it  is  convinient  to  set  W = col(E,  H),  and 


Q = 


c 

167T 


(5) 


where  qx  is  the  antisymmetric  tensor  dual  to  q (qxE  = q x E).  For  a time-harmonic  field,  the 
normal  component  Sq  of  the  time  average  Poynting  vector  S can  be  written  as  Sq  = q • S = 
W^QW.  Therefore,  the  condition  (Wn|<5|Wn)  = Nq  is  in  fact  the  normalization  to  the  beam 
energy  flux  Nq  through  the  plane  <tq: 

(Wn|0|Wn)=  f Sqd(TQ  — Nq. 


(6) 
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We  assume  here  that  the  tangential  component  t(6)  = k(6)  — q[q  • k(6)]  of  the  wave  vector 
k(6)  is  real  for  all  b E B,  and  that  the  mapping  b M-  t(fe)  is  one-one  (injective).  It  can  be  shown  [8], 
that  the  beams  W„  become  orthonormal,  if  B = Bu.  and  the  function  v{b)  is  given  by 


1 I NqJ{b) 

2tt  V g{b)Wi(b)QW{b)' 


(7) 


where  J(b ) = D{P)  / D(£l)  is  the  Jacobian  determinant  of  the  mapping  6 1(6),  calculated  in 
terms  of  the  local  coordinate  systems  (£*,  t = 1,2)  on  B and  (P,  j = 1,2)  on  the  t-plane, 
preserving  the  orientation  ( J(b ) > 0),  and  dB  = g{b)d^  d£2 . The  expression  under  the  square 
root  in  Eq.  (7)  has  to  be  finite  and  positive  almost  everywhere,  i.e.,  for  all  6 G B with  the 
allowable  exception  of  a set  of  measure  zero  in  B.  This  condition  can  be  met  with  appropriately 
chosen  beam  manifold  B,  operator  Q , and  normal  q.  In  some  cases,  when  B is  a proper  subset 
of  Bu,  i.e.,  B ^ Bu,  the  whole  set  of  beams  Wn  [Eq.  (3)]  can  not  be  orthonormalized;  instead, 
its  subset  may  be  orthonormalized. 

It  is  significant  that  the  beam  base  functions  K = K(6)  and  W = W(6),  the  beam  state 
function  v = z/(6),  and  the  orthonormal  functions  un  = un(b ) axe  defined  on  the  same  manifold 
Bu.  Therefore,  the  natural  coordinates  providing  the  coordinate  representation  of  un,  being  used 
as  integration  variables,  provide  also  the  natural  parametrization  of  the  beams  Wn(x)  [Eq.  (2)] 
or  Wn(r,  t)  [Eq.  (3)].  In  the  case  of  beams  defined  by  the  spherical  harmonics,  treated  in  Refs. 
[8,  9,  10],  Bu  is  a unit  sphere.  The  polar  angle  6 and  the  azimuthal  angle  <p  compose  the  natural 
coordinate  system  on  it.  These  coordinates  may  coincide  [8,  9]  or  be  closely  connected  [10] 
with  the  spherical  coordinates  of  the  propagation  direction,  but  this  is  by  no  means  always  the 
case.  In  particular,  to  compose  the  similar  beams  in  a biaxially  anisotropic  medium,  it  is  more 
advantageous  to  relate  9 and  ip  with  the  biaxial  coordinates  [11].  These  curvilinear  coordinates 
yield  the  parametric  representation  of  the  wavevector  surface,  which  makes  possible  to  obviate 
the  need  for  solving  algebraic  equations  to  describe  this  fourth  order  surface.  They  provide  also 
a very  convenient  means  for  decription  of  field  vectors  of  both  eigenwaves  and  beams  in  biaxially 
anisotropic  media. 

If  the  beam  Wn  (3)  consists  of  homogeneous  eigenwaves,  i.e.,  k*(6)  = k(6)  for  all  b € B,  it 
may  be  of  advantage  to  expand  it  into  a series  by  using  the  Rayleigh  formula 

+00  l 

e*kr  = 4tt  ]T  iljt(kr)  £ Yr'(k)Yr(r)  , (8) 

/— 0 m=—l 


where  k = k/A:,  r = r/r,  and  Y™  are  the  spherical  harmonics.  Substituting  the  expansion  (8) 
into  Eq.  (3),  we  obtain 


H-oo  l 

W(r, t)  = e+* £ y,m(f)W,m(r)  , (9) 

1=0  m=—l 

where 

W P(r)  = 4tt  [ jz(fe(6)r)yim*(k(6))l/(6)u(6)W(6)dB  . (10) 

J B 

Within  the  framework  of  this  description,  the  beam  is  characterized  by  a set  of  radial  vector 
functions  Wj71  = W;m(r).  In  an  isotropic  medium  (achiral  or  chiral),  these  relations  become 

4-oo  l 

W(r ,t)  = e-^£^(fcr)  £ 7^, 

l—Q  m=—l 


(ii) 
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where  the  coordinate  independent  vector  coefficients 

Wf1  = Att  [ y,m*(k(6))i/(6)«(^)W(6)dB  (12) 

‘ Jb 

completely  characterize  the  beam. 

4.  Conclusion 

New  families  of  exact  time-harmonic  solutions  of  Maxwell’s  equations  in  homogeneous  linear 
media  (isotropic,  chiral,  anisotropic,  and  bianisotropic)  and  free  space,  which  describe  orthonor- 
mal electromagnetic  beams,  are  obtained.  Owing  to  the  orthonormality  conditions,  these  beams 
form  convenient  functional  bases  for  more  complex  fields  and  provide  a helpful  technique  for 
modelling  the  beams  now  in  use  and  investigating  their  scattering  and  propagation  in  various 
media.  The  presented  results  form  a basis  for  generalizing  the  wave-splitting  technique  devel- 
oped in  [12, 13]  to  the  case  of  beams  in  plane-stratified  media.  They  can  be  used  for  modelling  of 
incident  and  scattered  beams  in  investigating  possible  physical  phenomena  caused  by  anomalous 
propagation  in  homogeneous  and  helicoidally  nonhomogeneous  bianisotropic  slabs. 

In  the  poster  presentation,  we  illustrate  the  obtained  solutions  by  calculated  spatial  distri- 
butions of  energy  density,  time  average  Poynting’s  vector,  and  local  polarization  parameters  for 
a number  of  beams  in  complex  media  and  free  space. 
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